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Abstract 

Performance of wave-energy devices of the oscillating water column (OWC) type is 
greatly enhanced when a resonant condition with the forcing waves is maintained. 
The natural frequency of such systems can in general be tuned to resonate with a 
given wave forcing frequency. In this paper we address the tuning of an OWC sea- 
water pump to polychromatic waves. We report results of wave tank experiments, 
which were conducted with a scale model of the pump. Also, a numerical solution for 
the pump equations, which were proven in previous work to successfully describe its 
behavior when driven by monochromatic waves, is tested with various polychromatic 
wave spectra. Results of the numerical model forced by the wave trains measured in 
the wave tank experiments are used to develop a tuning criterion for the sea-water 
pump. 
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1 Introduction 



Conversion of ocean wave energy into a useful form has been studied intensely 
duri ng the past two-decades, particu larly regarding generation of electricity 
(e.g. iRossl . Il995l ; IClement et all |2002| ) . Special attention has been given to one 
clas s of wave ene rgy converters based on an oscillating water column (OWC) 



see 



FalnesJ, I2OO2J, for a review of the OWC physics). In such systems, seawater 



is forced by incident waves to oscillate inside a partially-submerged chamber 
or duct. The free surface within this chamber compresses a volume of air which 
is, in most cas es, directed to flow t hrough a turbine t hus driving an electric 



generator (e.g. Falcao et al. . 1996 ; Heath et al. . 2000l ). Most R&D efforts 



on the subject h ave been focused on the wave energy absorption efficien cy 
optimisation (e.g. lEvansl . Il98ll ; iKordd . Il999l ; iPerdigao and Sarmentd . 120031 ). 



Air compression chamber 



Adjacent 
compartment 



Resonant duct 




Fig. 1. Schematic diagram of the seawater pump (from ICzitrom et all l2000hh 



An OWC wave energy device intended for sea -water pumping, invol ving no 
generation of electricity, has been described by ICzitrom et al. (boOOal lbh. The 
system has potential for various coastal management purposes, such as aqua- 
culture, flushing out of contaminated areas or the recovery of isolated coastal 
lagoons as breeding grounds. The sea-water pump (shown in figure 1) is com- 
posed of: a resonant duct, in which the wave oscillation, compared to the 
wave motion, is magnified and rectified by overtopping of the water column; 
an exhaust duct, through which water flow resulting from wave rectification 
is channeled to a receiving water body; and a variable volume air compression 
chamber coupling the two ducts. The wave-induced pressure at the mouth of 
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the resonant duct drives an oscillating flow that spills water into the com- 
pression chamber, and thence through the exhaust duct, with each passing 
wave. Air in the chamber behaves like a spring against which water in the 
resonant and exhaust ducts oscillates. Maximum efficiency is attained at res- 
onance when the system natural freque ncy of osc i llation coincides with the 



frequency of the driving waves (see e.g. iLighthilll . Il979l ). A resonant condi- 
tion can be obtained for different wave frequencies by means of a variable 
volume compression chamber that adjusts the hardness of the air spring. Re- 
sp qnse of the sys t em to monochromatic waves has been previously studied 



by ICzitrom et al.l (l2000bl ). They developed a tuning algorithm that predicts 
the optimal volume of air for the compression chamber, for given values of 
wave period and amplitude and tidal elevation. In this paper we address the 
problem of tuning the wave pump when forced by polychromatic waves. We 
report on wave tank experiments with a scale model of the pump as well as on 
the results of numerical simulations driven by real wave tank data. A tuning 
criterion is derived to optimize the system performance in these conditions. 

This paper is organized as follows: in section 2 we recall the equations describ- 
ing the seawater pump as well as previous results on the tuning to monochro- 
matic waves. Results of the experiments with a scale model and polychromatic 
forcing are reported in section 3. Tuning to polychromatic waves using exper- 
imental data and a numerical model of the pump is discussed in section 4 
followed by concluding remarks in section 5. 



2 Background 



2.1 System Equations 



The hydrodynamic theory for the OWC seawater pump studied here is based 



on what has been referred to as a heaving rigid-piston model (lEvansl . 11981 



Brendmo et al.l . Il996l ; iFalned . |2002| ). where the free surface inside the reso- 



nant duct is assumed to behave as a massless rigid piston. When the incoming 
waves are long with respect to the interior cross-section of the resonant duct, 
the deformation of the internal free-surface is negligible and the OWC is well 
described by the heave velocity of the 'rigid-piston'. The prediction of the sys- 



tem behavior is thus equivalent t o that of more general models (e.g. lEvans 



19821 ; ISarmento and Falcaol . Il985l ) that take into account the deformation of 



the int ernal free-surface. The pumping system was shown by ICzitrom et al 
( 12000al ) to be accurately described by equations (1-2) for the evolution of the 
water surface elevation inside each duct, when there is no spilling in the com- 
pression chamber. These equations were derived applying the time- dependent 
form of Bernoulli's equation to streamlines that join the free surfaces inside the 
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resonant and exhaust ducts wit h the surface s of the forcing and receiving wa- 
ter bodies, respectively (see e.g. iFalnesl . 12002 . cha pter 4). Terms were added to 
account for viscous loss es due to friction (se e also Perez et all Il996l ) and vor- 



tex formation (see also IKnott and Mackleyl . 119801 ). which are not considered 
in Be rnoulli's equation, and radiation damping (see also IKnott and Flowerl . 
1980). At the time of pumping, fluid surges from the resonant duct and spills 
into the exhaust side of the compression chamber. The fluid surface above the 
resonant duct bulges upward inducing a back pressure on the fluid in the duct, 
proportional to the bulge height. This effect has to be taken into consideratio n 



when solving the equations throughout the full cycle (jCzitrom et all l2000al ) 
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Equations (1-2) are written in non-dimensional form by scaling length and 
time with the amplitude A and frequency Q of an average forcing wave so 
that the forcing in equation (1) reads 



W(t) 



fj 



AA fl Q? 



sin (fit) . 



(5) 



A magnifying factor Af for the oscillation amplitude in each duct, with re- 
spect to the forcing wave, is considered. Subscripts 1, 2 and c in equations 
(1-2) correspond to the resonant and exhaust ducts and compression chamber 
respectively, x is the surface displacement, in either duct, relative to its' equi- 
librium position in the compression chamber, and L and A are lengths and 
areas respectively. L x and L 2 , the resonant and exhaust duct lengths, corre- 
spond to the physical len gth of the duc t s plus an added length due to edge 



effects at the mouths (see lCzitrom et alll2000bl ). In addition, P A is the atmo- 



spheric pressure, p is the sea-water density, g is the gravitational acceleration, 
7 = C p /C v = 1.4 is the ratio of the specific heat capacities of compressed air, 
Td is the height of sea level above the receiving body of water (including tides), 
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Vq is the compression chamber volume, H is the height of the equilibrium po- 
sition of the free surface in the compression chamber above the receiving body 
of water and C£>\,2 is a coefficient of nonlinear losses including vortex forma- 
tion, friction and radiation damping. Li and T<j are referred to the equilibrium 
position on the exhaust side of the compression chamber. 



The inertial term is first in both equations (1-2), where the water masses in 
the ducts are written in terms of the duct lengths — equation (3). The second 
and fifth terms come directly from Bernoulli's equation whereas the third term 
is added to represent the nonlinear losses due to vortex formation at the duct 
mouths, friction and radiation damping. Coupling of the two equations occurs 
through the air compression term (fourth in both equations) and the wave 
forcing is included on the right hand side of the resonant duc t equation (1 



A dime nsional analysis of the pump equations performed by ICzitrom et al 



(l2000al ) for a full-scale oceanic application, shows that the non-linear terms 
in equations (1-2) are relatively small, and that their main contribution is to 
limit flow through the pump. They showed that the pump behaves basically 
like a two- mass two-spring linear system (figure 2), where the two masses 
correspond to the water in the two ducts and the springs represent gravity and 
the air compression chamber. The restoring force on the oscillating system is 
thus provided by the compression of the air chamber and by the gravitational 
force. Non-linear losses, albeit small, are represented by the pistons. 
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Fig. 3. Modelled and observed resonant and exhaust duct water heig hts at the com- 



pressio n chamber and driving wave height plotted against time (from lCzitrom et al 
2000bh 



2.2 Numerical model 



A numerical model based on a finite-di fference integra t ion of the pump equa- 
tions (1-2), was thouroughly tested by ICzitrom et al.l (l2000bl ) for monochro- 
matic forcing waves. They compared output from numerical model runs for 
the free surface position in the ducts and flow through the pump, with experi- 
mental results obtained in wave tank experiments with a physical model of the 
pump. They used the pressure signal measured at the duct mouth as forcing 
instead of the sinusoidal wave of equation (5) and the system behaviour was 
very well predicted (see figure 3). 



2.3 Tuning to Monochromatic Waves 



An expression for the natural frequencies of oscillation of the sea-water pump 
when there is no spilling in the compression c hamber can be derived from the 
linearized version of equations (1-2). It reads (ICzitrom et all l2000bl ): 



I A ± VA 2 - 4B 



(6) 
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The two frequencies in equation (6) correspond to different modes of oscilla- 
tion of the water masses in the ducts. The higher frequency u>q, is associated to 
a synchronous oscillation of the two free surfaces that compresses and decom- 
presses the air chamber, and is the natural frequency excited for pumping. The 
second mode, with the lower frequency Uq , corresponds to a bodily oscillation 
of the two masses about their center of mass. 



Since expression (6) was derived from non-dimensional equations, the resonant 
condition useful for pumping occurs at u)q = 1, and the forcing wave frequency 
appears as Q in coefficients A and B. Solving for Vo provides the following 
simple analytica l algorithm for tuning the system for resonance at the given 



wave frequency (jCzitrom et all l2000al ) 
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This tuning algorithm, however, is not completely adequate when there is 
spilling in the compression chamber and the system is fully operational. Dur- 
ing pumping, the fluid that spills from the resonant duct ceases to weigh on 
the rest of the fluid in that duct, so that the gravitational spring restoring 
force becomes constant. This is equivalent to a softening of the combined air- 
cushion-gravity restoring force, since springiness is now only provided by the 
compression chamber. Pumping thus decreases the natural frequency of oscil- 
lation of the system so that, in order to maintain a resonant condition at the 
driving wave frequency, the compression chamber spring must be hardened, 
by decreasing its' volume. That is, resonance occurs at a volume of air smaller 
than the linear prediction Vquu when there is pumping. 



The numerical model was used by ICzitrom et al. fl2000bh to obtain a tuning 



algorithm for monochromatic waves for the fully operational system. This al- 
gorithm renders an optimal volume of air for the compression chamber (Vo), in 
terms of the wave period and amplitude and the tidal elevation. The resonant 
air chamber volume (Vo) was thus found to increase with increasing wave pe- 
riod and decreasing tidal height. They also showed that the flow rate through 
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Fig. 4. Schematic diagram of the scale model tested in the wave tank (see text). 

the pump at resonance (Q ) increases with increasing wave period and wave 
amplitude. 



3 Experiments 



In this section we describe wave tank experiments conducted with a scale 
model of the pump driven by polychromatic waves. The model was designed 
on a 1:20 scale with respect to a full-scale oceanic application a nd is similar to 



that us ed for the monochromatic wave experiments reported in ICzitrom et al. 



( I2000bl ). A schematic diagram of the experimental setup is shown in figure 4. 



The mouth of the resonant duct (A) was placed in the wave tank exposed to 
the passing waves, while the exhaust duct (B) was placed in a bucket of water 
(C), which was leveled flush with the surface of the tank. Water spilt from the 
bucket was collected in a recipient (D). The pump was primed at the start of 
the experiments by creating a partial vacuum that brought water up from the 
tank and bucket, to a working level in the compression chamber (E). We define 
the portion of resonant duct extending above this level as the sill height. The 
volume of air in the compression chamber was modified by changing the level 
of water in the connecting chamber (F); through interchange with the storage 
tank (G). The model was fully instrumented with water height sensors of the 
capacitance type in the resonant and exhaust ducts, as well as in the water 
deposit (D), used to monitor flow through the pump continuously. A water 
height sensor was also placed at the wave tank surface, above the resonant 
duct mouth, to measure the driving wave height. The sensors were sampled at 
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0.1s intervals. Dimensions of the various components are explicited in Table I. 



u 


D 1 


A l 


U 


D 2 


A 2 


5.25m 


0.056m 


0.00985m 2 


16m 


0.036m 


0.00407m 2 



D c 


A c 


V 


Sill height 


0.145m 


0.06605m 2 


0.033m 3 -0.073m 3 


0.005m-0.03m 



Wave tank depth 


Resonant duct mouth depth 


0.86m 


0.279m 



Table I. Dimensions of the various components of the model pump. 

The model pump was driven with waves of various spectra, which were de- 
signed to mimic real sea surface waves and also to allow for differentiation of 
the response of the system to each frequency component. For the data shown 
in this paper, the wave generator was programmed to include four different 
wave periods (1.8, 2, 2.2 and 2.4 seconds) while keeping the amplitude at 0.04 
m for all components. Traces in time obtained from the different sensors for 
one of the experiments are shown in figure 5. 

It can be seen in this figure that the pulsating forcing wave drives the response 
in the resonant duct and in the flow. The regions with largest slope in the flow 
trace lag behind the largest oscillations in the resonant duct which in turn lag 
behind the driving wave pulsations. From this phenomenon we can charac- 
terize a response time for the experimental system, between the largest wave 
oscillations and the largest observed flow rate, of approximately 12 seconds. 
It can also be seen that the events of large oscillations in the resonant duct, 
which produce pumping, are associated initially with an increased water level 
in the exhaust duct side of the compression chamber, which subsequently os- 
cillates with a period greater than the driving wave. This lower frequency is 
related to the second natural frequency of oscillation of the sea-water pump, 
identified as Uq in expression (6), and is excited by the change in level which 
results from spilling in the compression chamber. 

In figure 6, Fast Fourier Transforms (FFT) of 1024 data points of the forcing 
wave signal and the oscillations within the ducts are shown. Comparing figures 
6. a and 6.b, it can be seen that, for the resonant duct, the main driving wave 
frequencies are transformed so that the resonant frequency peak is clearly 
larger than the others for that particular experiment. The lower frequencies 
also present in the resonant duct response correspond, on the one hand, to 
the second natural frequency of oscillation noted in the previous paragraph, 
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Fig. 5. Time series of the water height measuremments in the wave tank, the reso- 
nant and exhaust ducts (in metres) and the flow recipient (in cubic metres). 

and to combination frequencies nfi — mfj where n and m are integers (n = 
m = 1 most important) and and fj are any of the main driving frequencies. 
These combinations are to be expected considering that the frequency of the 
pulsations in the driving wave must correspond to differences between the 
main forcing frequencies. In the exhaust duct (figure 6.c), these low frequency 
components are the most important part of the oscillation. 



4 Tuning to Polychromatic Waves 



4-1 FFT analysis 



In figure 7, we show FFT's of the forcing wave signal and of the oscillations 
in the two ducts measured in four different cases. These four experiments 
were conducted in the same wave conditions (a polychromatic wave with four 
components T po i y = 1.8,2,2.2 and 2.4) and only the volume of air in the 
compression chamber had a different value in each case, corresponding to the 
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Fig. 6. Fast Fourier Transforms of (a) the driving wave height, water height in 
the (b) resonant duct and (c) exhaust side of the compression chamber. These 
FFTs correspond to the time series shown in figure 5 and are all normalized to the 
maximum of (a). 



prediction Voun of equation 10 for each period component. Also shown in 
these figures are the two natural frequencies of oscillation (f + and f _ = 1/T) 
obtained from equation (6) for each volume of air in the compression chamber. 
The energy in the forcing wave (O's in figure 7), which is injected mainly in 
the four principal peaks of the signal, is redistributed in the response of the 
system. The frequency content in the resonant duct signal (D's in figure 7) is 
clearly different for each case and it can be seen that the largest peak occurs 
near the natural frequency of oscillation f + predicted by the linearized model. 
It is apparent that the system responds with greater strength to the frequency 
nearest its natural frequency of oscillation, despite the polychromatic nature 
of the forcing wave. 
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Fig. 7. Fast Fourier Transforms of experimental data of the forcing wave (O) and 
the oscillations in the ducts (□: resonant duct signal; +: exhaust duct signal, very 
weak with respect to the driving wave and the resonant duct). The same driving 
wave of 4 period components T po i y = 1.8,2,2.2 and 2.4 was used in all cases. The 
volume of air in the compression chamber Vq was varied in the four cases shown to 
match the value determined by each value of VoLm(^poiy) ( see text). Its value in m 3 
was (a) 0.0058, (b) 0.0074, (c) 0.0091 and (d) 0.011. The two natural frequencies of 
oscillation for each configuration are shown as vertical lines. 



4-2 Tuning from the experimental data 



The air chamber volume for which flow through the pump is greatest for each 
particular configuration (henceforth referred to as the resonant volume), can 
be determined from the experimental flow data. Plots of average flow rate 
as a function of the volume of air in the compression chamber for various 
sill heights are shown in figure 8, for the same 4-component forcing wave of 
figure 7. The long vertical lines show the linearized model prediction of the 
resonant air volume {Voun) for each forcing wave frequency component (the 
prediction for the 1.8 seconds component lies below the range covered in this 
plot). The vertical line furthest to the right corresponds to the largest wave 
period present (2.4 s). If we consider the monochromatic responses, according 
to the algorithm for resonant flow (Qo) mentioned above and given that all 
four main components had essentially the same amplitudes, we should expect 
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Fig. 8. Average flow rate obtained in the experiments vs. volume of air in the 
compression chamber for the same forcing wave of figure 7. The long vertical lines 
show the linearized model prediction Von n for each component. The short vertical 
lines show the numerical model prediction Vq for the greater period and each sill 
height (given in meters). 

the response to the component of the largest period to produce the greatest 
pumped flow, fn figure 8, the observed maximum flow rates occur for volumes 
of air in the compression chamber bounded by the linearized model predic- 
tions (VoLin) for the two largest periods. This suggests that the components 
of larger periods determine the collective resonant response of the system to 
polychromatic waves. It can also be seen in figure 8, that the resonant vol- 
ume, as denned by the maximum average flow rate is, for all sill heights, always 
smaller than the Voun for the maximum period component. Th is is consistent 



with the monochromatic experiments of ICzitrom et al.l (l2000bl ) who observed 
that pumping shifts the actual resonant volume to smaller values than the lin- 
earized model prediction. Similar to the monochromatic case, greater pump- 
ing induces a greater shift towards smaller resonant air chamber volumes. The 
short vertical lines in figure 8 are the values for the resonant volume obtained 
from monochromatic numerical model experiments with the greatest period 
and six sill heights. These were obtained by running the model for a series of 
compression chamber air volumes covering the range shown in figure 8. It is 
clear that the resonant air chamber volume predicted by the monochromatic 
model closely resembles the observed volume that produces maximum for each 



13 



1.2 



1 



0.8 



0.6 




0.4 



0.2 











0.2 



0.4 



0.6 



0.8 



1 



Frecfjency(Hz) — WiveFFT — Op 



Fig. 9. Fast Fourier Transform of the forcing wave signal and the corresponding 
graphic for Qo(T) for a two-component forcing wave (of 2 and 2.6 seconds periods). 
Each curve is normalized with respect to its maximum. 

case. This result encourages us to build a tuning criterion for waves with vari- 
ous frequency components based on the numerical monochromatic algorithms 
for V and Q . 

It should be noted that a forcing wave with the same amplitude for all of its' 
components is, of course, unlikely a common real sea-surface wave. Thus, the 
greatest pumping can occur when the system is tuned to a component that is 
not necessarily the one with the largest period. In the next section we use the 
flow prediction obtained numerically for monochromatic waves (Qo), which 
takes into account the wave amplitude for a given frequency, and use it to 
choose an optimal tuning condition. 



4-3 Numerical experiments with polychromatic waves 

We performed a series of numerical simulations using the same configuration 
of the wave tank experiments, but driving the model with different synthetic 
waves of multiple frequency components. The mathematical expression of the 
forcing waves used in the numerical experiments was: 



A i sin 



2nt 
~~T~ 
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Fig. 10. Flow vs. volume of air in the compression chamber obtained with the 
numerical model forced by: a signal with components of 2 and 2.6 s (Q(Ti,T2)); 
only the 2s component (Q(Ti)); only the 2.6s component (<2(?2)). 

We report on two particular cases that let us discuss the basic features of 
the polychromatic response: In the first case, the forcing signal was a two- 
component wave (of Ti = 2 and T 2 = 2.6 seconds) whose FFT is shown in a 
heavy line in figure 9. We used a greater amplitude for the higher frequency 
deliberately A 2 = 0.6A l5 to show that the wave amplitude must also be used 
in determining the maximum flow through the pump. The light line in the 
same figure is the monochromatic flow prediction Q , calculated with the 
frequency / = T _1 and corresponding amplitude values from the FFT. It 
can be seen that, in contrast with the experimental observations which were 
conducted with the same amplitude for all frequency components, here we 
find that the smaller period produces a greater flow. Clearly, the amplitude 
of each component must be taken into account to choose the most efficient 
tuning condition. We remark, however, that the ratio between the expected 
maximum flows (Q in figure 9) for the two frequencies, is noticeably closer to 
1 than the ratio between the amplitudes of each component in the FFT of the 
forcing wave. This shows that energy from components with larger periods is 
used more efficiently by the system. 

Once we have chosen the wave frequency component which produces the great- 
est flow, the appropriate air chamber volume can be calculated using the 
monochromatic numerical algorithm for Vq. In order to test the effectiveness 
of this criterion as a selection tool for the optimal volume of air of the compres- 
sion chamber, in figure 10 we show the flow rate through the pump plotted 
against volume of air in the compression chamber, computed with the nu- 
merical model forced by the wave signal shown in figure 9 (Q(T 1 ,T 2 )). Also 
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Fig. 11. Fast Fourier Transform of the forcing wave signal and the corresponding 
graphic for Qo(T) for a four-component wave (of 2.9, 2.75 and 2.6 and 2 seconds). 
Each curve is normalized with respect to its maximum. 

shown in this figure (Q(Ti) and Q(T 2 )) are the flow rates obtained by forcing 
the numerical model with monochromatic waves of frequencies and amplitudes 
corresponding to each one of the components in figure 9. The abscissa in figure 
10 is normalized with the optimal volume of air predicted by Vq for the higher 
frequency. The flow response for the polychromatic forcing (Q(T 1 ,T 2 )) has 
two maximums, which correspond to the resonant volumes for each frequency 
component. It is confirmed that the absolute maximum corresponds to the one 
predicted by Qq in figure 9 (1 in the horizontal scale). The correspondence of 
the maximum flows in the poly and monochromatic wave forcings in figure 10 
shows that the attainable flow with one frequency is essentially not affected 
by the presence of the other. That is, for this separation between frequencies, 
the pump behaves similar to a flute which amplifies each frequency distinctly, 
depending on the natural frequency of oscillation to which it is tuned. 

We now analyze another case in which the forcing wave has four components 
of periods Tj = 2,2.6,2.75 and 2.9 seconds. As it can be seen in figure 11, 
the signal is composed of a block of three frequencies close to each other (cor- 
responding to the wave periods of 2.9, 2.75 and 2.6 seconds) plus a fourth 
frequency, slightly more separated, corresponding to 7\ = 2s. As with the 
previous case, we attempt using the Q curve to determine which compo- 
nent we should choose for tuning the system. If we guide ourselves only by 
the maximum flow value, we would again choose the 2 seconds component 
that is separated from the block. Similar to figure 10, we used the numerical 
model to test the validity of this tuning criterion as shown in figure 12. Here 
it can be seen that, while the absolute maximum flow rate in the polychro- 
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Fig. 12. Flow vs. volume of air in the compression chamber obtained with the numer- 
ical model forced by: a signal with components of Tj = 2, 2.6, 2.75 and 2.9 seconds 
(solid line, Q(Ti, . . . T 4 )) and four series obtained with monochromatic forcings for 
each component (Q(Tj)). 



matic response (Q(T±, . . . T 4 )) is correctly identified with the corresponding 
monochromatic response (Q(Ti)), the flow values Q(Ti, . . . T 4 ) for the air vol- 
umes resonant with the lower frequencies are clearly higher than the maxima 
of the monochromatic responses (Q(T 2 ), Q(T 3 ) and Q(T 4 )). This shows that 
the responses to the three lower frequencies in the forcing wave interact to give 
a larger flow than expected from each component separately and indicates the 
existence of a bandwidth of re sonance. The bandwidth in oscillating systems 
is related to damping (see e.g. iFalned . 120021 ) . For the seawater pump damping 
is dominated by the pumping process itself : when the amplitude of oscilla- 
tion increases in the resonant duct, water is spilt at the sill during a larger 
part of the period. This determines that, in the pumping regime, the limit im- 
posed on the amplitude by the pumping process dominates over other damping 
mechanisms such as friction or vortex formation. The resonance bandwidth of 
the pumping system can be denned using different norms. We choose it to 
be the period interval for which flow through the pump is at least 80% its 
maximum value, and estimate it using flow data (not shown here) from the 
scale model experiments with different monochromatic forcing periods. For 
the experimental configuration of figure 12, a resonance bandwidth covering 
a broad period interval between 1.5 and 3 seconds is obtained. This rather 
large bandwidth can explain the previously observed fact that the polychro- 
matic response Q(Ti, . . . T 4 ) in figure 12 shows greater pumped flow, at the air 
volumes resonant with the lower frequencies, than any of the monochromatic 
responses (Q(T 2 ), Q(T 3 ) and Q(T 4 )). 
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Fig. 13. Filtered FFT of the forcing wave signal of figure 11 with a running mean of 
0.04 Hz and the corresponding Qo(T) curve; each curve is normalized with respect 
to its maximum. 



In real sea-surface waves, some or all of the frequencies present in a block are 
more likely to persist in time, than a single frequency on its' own. Furthermore, 
because of the resonance bandwidth discussed in the previous paragraph, a 
better pump performance can be expected when exploiting the collective re- 
sponse of the system to various peaks of a block. An ultimate tuning criterion 
for irregular sea waves should therefore take into account information on the 
frequency-domain distribution of the forcing. We propose using a running 
mean of a given FFT as a method of pondering the relative importance of 
blocks of frequencies versus single frequencies. Such a running mean would 
diminish the amplitude of single peaks, while allowing blocks to better retain 
their relative weight. As an example of how this procedure could work, in 
figure 13 a running mean of the FFT in figure 11 with a width of 0.04 Hz is 
shown (we choose the width for the filtering as 10% of the estimated band- 
width, which was ps 0.4Hz for this case). As expected, the reduction of the 
single frequency peak is noticeably greater than that of the three frequency 
block. Also shown in this figure is the expected flow Qo as calculated using 
the filtered FFT amplitudes. It is clear that a somewhat greater flow is now 
associated to the block of frequencies than to the single frequency peak. 

Using this type of filtering of the frequency content of the forcing wave signal 
can be an important tool in the design of the tuning algorithm for the sea- 
water pump because, as mentioned before, in a constantly changing sea, blocks 
of frequencies in the forcing spectrum have the added convenience of a likely 
greater persistence in time than single frequency peaks. 
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5 Conclusions 



The response of an oscillating-water-column (OWC) seawater pump to poly- 
chromatic waves was studied using wave tank tests with a scale model and 
numerical simulations. We show that, in spite of being driven by polychromatic 
waves, the system responds mainly in a monochromatic way, amplifying a se- 
lected frequency from the wide spectrum present in its forcing. The measured 
flow through the pump as a function of the experimental tuning parame- 
ter (the volume of air in the compression chamb er) retrieves th u s the s ame 
basic qualitative features that were reported by ICzitrom et al.l (J2000bJ) for 
experiments with monochromatic waves. The main observations are that the 
natural frequency of oscillation of the OWC decreases due to the overtopping 
associated with the pumping process, and that greater pumping is obtained 
when tuning the system to wave components of lower frequencies, given equal 
amplitudes. For the experiments reported here, the observed resonant volume, 
defined as the volume of air in the compression chamber for which flow through 
the pump is greatest, is predicted remarkably accurately by monochromatic 
algorithms using the lowest frequency peak present in the forcing spectrum. 
Numerical experiments with synthetic forcing waves showed nonetheless that 
the polychromatic response of the system is not always merely a superposition 
of monochromatic responses. A cooperative interaction between the responses 
to slightly different frequencies present in the forcing wave was evidenced, 
confirming the existence of a resonance bandwidth and hinting to the use 
of a frequency-space-averaged signal as input of the proposed polychromatic 
tuning algorithm. The latter is based on the monochromatic algorithms for 
resonant volume of air in the compression chamber (Vo) and resonant flow rate 

(Qo). 



The tuning process proposed for the sea-water pump in an irregular sea con- 
sists of various steps that can be easily performed by a pre-programmed mi- 
crochip and is summarized in figure 14. First, the FFT analysis would be 
performed on a series of measurements taken with a wave sensor, and the 
resulting spectrum would be filtered to ponder the frequency content of the 
signal. The maximum expected flow would then be identified using the Qq 
calculation from the filtered amplitudes and corresponding frequencies and, 
finally, the resonant volume of air for optimal sea-water pump performance 
would be computed using the Vo algorithm with the appropriate frequency 
and unfiltered amplitude of the chosen component. This tuning procedure 
thus outputs the volume of air in the compression chamber that will give the 
best pump performance. An actual control system modifying the volume of air 
in the compression chamber could then perform its task at any given time tak- 
ing only as input a series of measurements of the wave signal at the resonant 
duct mouth. 
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Fig. 14. Tuning protocol diagram (see text for details). 
Acknowledgments 



We wish to thank Tony Dalrymple, Joe Hammack, Nobu Kobayashi and Brad 
Johnson for allowing and helping us to use the wave tanks at the Center for 
Applied Coastal Research, The University of Delaware. Thanks are also due 
to Ranulfo Rodriguez Sobreyra for developing and implementing the sensing 
and data acquisition system used during the experimental work, and Esteban 
Prado, Catalina Stern and Arturo Olvera for useful discussions. We greatly 
acknowledge the John D. and Catherine T. MacArthur Foundation for their 
generous support through a grant from the Fund for Leadership Develop- 
ment. Support was also provided by the Direction General de Apoyo al Per- 
sonal Academico (Projects IN-104193, 106694, 107197 and a Graduate Stud- 
ies Scholarship for RGD) and the Instituto de Ciencias del Mar y Limnologfa 
(Project "Energfa de Oleaje" No. 139), both from the National University of 
Mexico (UNAM). Support was also provided by CONACyT through project 
'Matematicas no-lineales en la fisica y en la ingenieria' (FENOMEC), No. 
625427-E. 



References 

Brendmo, A., Falnes, J., Lillebekken, P., 1996. Linear modelling of oscillating 
water columns including viscous loss. Appl. Ocean Res. 18, 65-75. 

Clement, A., McCullen, P., Falcao, A., Fiorentino, A., Gardner, F., Ham- 
marlund, K., Lemonis, G., Lewis, T., Nielsen, K., Petroncini, S., Pontes, 
M., Schild, P., Sjostrom, B., S0 rensen, H., Thorpe, T., 2002. Wave energy 
in europe: current status and perspectives. Renewable Sustainable Energy 
Rev. 6, 405-431. 

Czitrom, S., Godoy, R., Prado, E., Perez, P., Peralta-Fabi, R., 2000a. Hydro- 



20 



dynamics of an oscillating water column seawater pump, part I: Theoretical 
aspects. Ocean Eng. 27, 1181-1198. 

Czitrom, S., Godoy, R., Prado, E., Olvera, A., Stern, C, 2000b. Hydrody- 
namics of an oscillating water column seawater pump, part II: Tuning to 
monochromatic waves. Ocean Eng. 27, 1199-1219. 

Evans, D., 1981. Power from water waves. Annu. Rev. Fluid Mech. 13, 157- 
187. 

Evans, D., 1982. Wave power absorption by systems of oscillating surface 
pressure distributions. J. Fluid Mech. 114, 481-499. 

Falcao, A., Sabino, M., Whittaker, T., Lewis, A., 1996. Design of a shoreline 
wave power pilot plant for the island of pico, azores. In: Proc. 2nd European 
Wave Power Conference. Lisbon, Portugal, November 1995, pp. 87-93. 

Falnes, J., 2002. Ocean waves and oscillating systems. Cambridge University 
Press, Cambridge, UK. 

Heath, T., Whittaker, T., Boake, C, 2000. The design, construction and op- 
eration of the limpet wave energy converter (islay Scotland). In: Proc. 4nd 
European Wave Energy Conference. Aalborg, Denmark, pp. 49-55. 

Knott, C, Flower, J., 1980. Measurement of energy losses in oscillatory flow 
through a pipe exit. Appl. Ocean Res. 2 (4), 155-164. 

Knott, C, Mackley, M., 1980. Eddy-motions near plates and ducts induced 
by water waves and periodic flows. Phil. Trans. Royal Soc. 294A (1412), 
599-623. 

Korde, U., 1999. On providing a reaction for efficient wave energy absorption 
by floating devices. Appl. Ocean Res. 21, 235-248. 

Lighthill, J., 1979. Two-dimensional analyses related to wave-energy extrac- 
tion by submerged resonant ducts. J. Fluid Mech. 91, 253-317. 

Perdigao, J., Sarmento, A., 2003. Overall-efficiency optimisation in owe de- 
vices. Appl. Ocean Res. 25, 157-166. 

Perez, P., Czitrom, S., Peralta-Fabi, R., Treviho, C, 1996. Friction in os- 
cillating flows. In: Proc. 2nd European Wave Power Conference. Lisbon, 
Portugal, pp. 375-380. 

Ross, D., 1995. Power from the waves. Oxford University Press, Oxford, UK. 

Sarmento, A., Falcao, A., 1985. Wave generation by an oscillating surface 
pressure and its application in wave-energy extraction. J. Fluid Mech. 150, 
467-485. 



21 



